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Abstract Using the water isotope- and vapor source distribution (VSD) tracer-enabled Goddard Institute for 
Space Studies ModelE-R, we examine changing El Nino-Southern Oscillation (ENSO)-like expressions in the 
hydrological cycle in a suite of model experiments. We apply strong surface temperature anomalies associated 
with composite observed El Nino and La Nina events as surface boundary conditions to preindustrial and 
mid-Holocene model experiments in order to investigate ENSO-like expressions in the hydrological cycle under 
varying boundary conditions. We find distinct simulated hydrological anomalies associated with El Nino-like 
("ENSO warm ") and La Nina-like ("ENSO C ool") conditions, and the region-specific VSD tracers show hydrological 
differences across the Pacific basin between El Nino-like and La Nina-like events. The application of ENSO C ool 
forcings does not produce climatological anomalies that represent the equal but opposite impacts of the 
ENSOwarm experiment, as the isotopic anomalies associated with ENSOwarm conditions are generally stronger 
than with ENSO C ool and the spatial patterns of change distinct. Also, when the same ENSO-like surface temperature 
anomalies are imposed on the mid-Holocene, the hydrological response is muted, relative to the preindustrial. 
Mid-Holocene changes in moisture sources to the analyzed regions across the Pacific reveal potentially complex 
relationships between ENSO-like conditions and boundary conditions. Given the complex impacts of ENSO-like 
conditions on various aspects of the hydrological cycle, we suggest that proxy record insights into paleo-ENSO 
variability are most likely to be robust when synthesized from a network of many spatially diverse archives, which can 
account for the potential nonstationarity of ENSO teleconnections under different boundary conditions. 


1. Introduction 

The El Nino-Southern Oscillation (ENSO) is characterized by sea surface temperature (SST) anomalies in the 
tropical Pacific that influence remote rainfall through atmospheric teleconnections and dominate climatic 
variability in many regions. The strength of remote ENSO teleconnections (i.e., the relative influence of a 
remote climatic mechanism on regional and local climates) is dependent on changes in the large-scale state 
of the coupled ocean-atmosphere system [Gallant et al., 2013], and potential changes in the amplitude, 
length, and frequency of ENSO events under future-changing boundary conditions remain poorly 
constrained [Collins etal., 2010]. In this present study, we investigate expressions of ENSO-like climatic modes 
in the hydrological cycle using a global climate model (GCM). As substantially different ENSO expressions 
may exist in the hydrological cycle under different orbital forcings, we examine the relationship between 
boundary conditions and changes in the hydrological cycle under ENSO-like forcings using the preindustrial 
and the mid-Holocene as research targets that characterize distinct background climatic states. 

1.1. Proxy S 1s O Reconstructions 

Paleoclimatic studies can provide a long-term context for understanding varying ENSO expressions [Jones 
and Mann, 2004; Gergis et al., 2006], although there remain obstacles to increasing our understanding of 
ENSO variability through both paleoarchives and model-based approaches. First, precisely dated, well- 
resolved, readily interpretable proxy records are limited. While oxygen isotope (5 1s O, 8 in permil units, %o, of 
the subscripted value relative to a known standard) archives can provide high temporal resolution records 
from ENSO-impacted areas [Cobb et al., 2003, 2013], proxy 6 ls O is often dependent on multiple, complex 
climatic S 18 0. Coral 8 1s O responds to water temperature changes together with rainfall variability, due to the 
influence of isotopically depleted rainfall on surface ocean composition [Tudhope et al., 2001]. Speleothem 
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and ice core-based § 18 0 records provide complementary climatic archives [ Lachniet et al., 2004], but these 
may integrate changes in temperature or precipitation, together with vapor source due to circulation 
changes [ Vuille et al., 2003]. 

In addition, the nonstationary nature of ENSO teleconnections can limit the utility of proxy-derived ENSO 
reconstructions [Gallant eta!., 2013]. As the spatial influence of ENSO likely changes under different boundary 
conditions, preserved climatic signals may become convoluted [Wilson et al., 2010] and interpolating climatic 
conditions from local reconstructions to basin-wide generalizations may be spurious [Rosenthal and Broccoli, 
2004], Hence, multiproxy data from many widely spaced regional proxies are most useful for obtaining 
precise information of temporal changes in large-scale ocean-atmosphere processes [e.g., Emile-Geay et al., 
2012], Overall, ENSO-sensitive proxy records indicate distinct changes in the nature of teleconnection 
patterns, together with differences in ENSO amplitude, frequency, and length compared to the present day 
[e.g., Tudhope et al., 2001 ; Bush, 2007; McGregor et al., 201 3a]. 

1.2. Previous Modeling Work 

Global climate models provide an additional tool for understanding ENSO variability on longer timescales 
and also for understanding how 8 1s O archives record climatic variability [ Toniazzo , 2006; Brown et al., 2008; 
Chiang et al., 2009; Conroy et al., 2013]. As both the climatic mean state and dynamics of the coupled ocean- 
atmosphere ENSO system in the tropical Pacific were likely different in the mid-Holocene compared with 
present day, this provides a useful period to examine changes in the ENSO system [Luan eta!., 2012], Previous 
model-based studies indicate that ENSO is sensitive to the orbital variations, and the dynamical interaction 
between the tropical atmosphere and oceans is particularly influenced by precessional modulation of the 
seasonal cycle [Clement et al., 1 999]. A reduction in ENSO amplitude in the mid-Holocene arises from changes 
in seasonality that influences the depth and intensity of the thermocline and the timing of the development 
of ENSO associated [Luan et al., 2012]. 

Modeling analyses that explicitly incorporate water isotope tracers also move us closer to making 
comparisons between proxy records and model data [Brown et al., 2006, 2008; Schmidt et al., 2007]. Water 
isotope tracer analyses have helped constrain 8 ls O anomalies associated with modern ENSO conditions and 
determined regions where simulated S 18 O p is significantly correlated with ENSO variability [Vuille and Werner, 
2005; Tindall et al., 2009], Model-based approaches help elucidate the relationship between ENSO and 
hydrologic variability, but they have not fully explored the characteristics of changes in ENSO expressions, 
including isotopic variability, under differing mean climate states. 

In this present study, we investigate modeled hydrological expressions of ENSO-like climatic modes around 
the tropical Pacific using multiple simulations of past and present climates. By employing the isotope- 
enabled Goddard Institute for Space Studies (GISS) ModelE GCM, we investigate the pattern of isotopic 
anomalies associated with forced ENSO-like events under varying boundary conditions. In addition, we 
incorporate a novel set of vapor source distribution (VSD) tracers as an additional diagnostic tool to assess 
changes in atmospheric circulation at proxy sites in present-day ENSO-sensitive regions. We use modeled 
mid-Holocene time slices forced by ENSO-like sea surface temperature anomalies and examine whether 
ENSO-associated hydrological expressions in the mid-Holocene necessarily follow those simulated during the 
preindustrial. Overall, this study provides a model-based understanding of the potential nonstationarity of 
ENSO-like hydrological expressions during the mid-Holocene compared to the preindustrial. 

2. Methods 

2.1. Model Description 

Simulations were made using the coupled atmosphere-ocean GISS (Goddard Institute for Space Studies) 
ModelE-R, which is the version submitted to the Coupled Model Intercomparison Project Phase 3 (CMIP3), 
with the addition of water isotope and tracers. The horizontal resolution is 4° latitude x 5° longitude with 20 
vertical levels up to 0.1 hPa in the atmosphere [Schmidt et al., 2006] and a 1 3 layer Russell ocean model of the 
same horizontal resolution [Hansen et al., 2007]. Water isotope tracers ( 1 H 2 16 0, "normal" water; 2 H 1 H 1s O or 
HDO, reported as 8D; and 'H^O, S 18 0) are incorporated into the atmosphere, land surface, sea ice, and 
ocean. Water isotopes are tracked through all stages of the hydrologic cycle and are advected throughout the 
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Table 1. Summary of Model Boundary Conditions for the Ok and 6k Experiments Analyzed Here, Including Greenhouse Gases [Indermuhle et at., 1999] and 
Perihelion in Julian Days [Berger and Loutre, 1 991 ] a 


Time (kyr) 

ENSO Phase 

C0 2 (Fraction) 

CH 4 (Fraction) 

N 2 0 (Fraction) 

Perihelion 

Prescribed Surface Temperatures 

0k 

NEUTRAL 

1 

1 

1 

2.85 

N/A 


WARM 

1 

1 

1 

2.85 

El Nino-like average surface temperature anomalies 


COOL 

1 

1 

1 

2.85 

La Nina-like average surface temperature anomalies 

6k 

NEUTRAL 

0.95 

0.71 

0.82 

263.42 

N/A 


WARM 

0.95 

0.71 

0.82 

263.42 

El Nino-like average surface temperature anomalies 


COOL 

0.95 

0.71 

0.82 

263.42 

La Nina-like average surface temperature anomalies 


a For each time slice, parallel experiments are also conducted with positive and negative ENSO-like surface temperature phases (ENSO W arm and ENS0 CO ol)- 
Sea surface temperatures for ENSO-like phases are derived from observational temperatures [ Reynolds and Smith, 1994] and NIN03.4 indices over the period 
of 1981 to 2003. 


model like water. At each phase change a fractionation is applied that explicitly determines equilibrium 
fractionation, with parameterizations accounting for kinetic fractionations [ Schmidt et at., 2005]. 

In addition to the water isotope module, we include a suite of water vapor source distribution tracers in the 
model so that sources of water vapor can be traced back through any cloud processes to the site of surface 
evaporation. The use of VSD tracers provides an additional circulation diagnostic for identifying regional- 
scale hydrological changes. The tracers are altered by atmospheric transport and condensation processes 
analogously to a nonfractionating water isotope tracer [Kelley, 2003]. The VSD is the integrated mass of water 
vapor in each model cell and is expressed as an area integral of evaporative input unique to that cell. This 
study uses spherical harmonics as VSD basis functions, which are not anchored to any particular geographic 
boundary and require no prior definition of regions. The precipitation source distribution is a subset of the 
VSD, defined where vapor condenses to liquid. We include 144 tracers that are resolved to wave number 11, 
providing an effective horizontal resolution to approximately 8° latitude x 10° longitude. It should be noted 
that factors such as land-sea contrasts cause real-world precipitation source distributions to vary unevenly 
over planetary scales, and hence, spatial VSD patterns cannot be interpreted at scales finer than their 
effective resolution. 

2.2. Experimental Design 

We conducted the experiments for this study in several stages (Table 1). As the VSD tracers are expensive to 
compute, an initial series of coupled atmosphere-ocean simulations were run without the VSD tracers 
[LeCrande and Schmidt, 2009]. These coupled simulations were run to equilibrium (100 s to 1000 years) for 
preindustrial (0 k) and mid-Holocene (6 kyr BP, "6 k") [LeGrande and Schmidt, 2009] mean climate conditions. 
The preindustrial control simulation (0 k) was run with all boundary conditions and atmospheric composition 
appropriate to the preindustrial period (circa 1880). In the 6 k experiment, greenhouse gas concentrations 
were adjusted according to ice core reconstructions [Indermuhle et a!., 1 999; Brook et al., 2000; Sowers et al., 
2003]; seasonal insolation was also changed relative to orbital forcings [Berger and Loutre, 1991], 

The second phase of our experiment includes the VSD tracers, but in atmosphere-only simulations where 
surface conditions (SST and sea ice) are determined from the described coupled simulations. The 
atmosphere-only, VSD-enabled model was run for 6 years, with the average of the final 5 years' results 
reported here. Using this approach, once the sea surface conditions are determined (from long runs to 
equilibrium using the coupled model), the atmospheric model can be driven for much shorter runs using 
the surface conditions determined from the coupled model. A preindustrial simulation with VSD tracers 
included was also conducted in fully coupled mode in order to test the validity of using results from 
atmosphere-only simulations, and this shows only small differences in precipitation source distributions, 
indicating that the atmosphere-only approach is valid. 

Next, we performed a set of ENSO-like simulations for each of our mid-Holocene and preindustrial 
experiments, in which either strong El Nino (warm) or La Nina (cool) surface temperature anomalies are 
imposed in the model in order to force an atmospheric response. We create a global composite of El Nino and 
La Nina anomalies using observational temperatures [Reynolds and Smith, 1994] and NIN03.4 indices over the 
period of 1 981 to 2003. We regress NIN03.4 indices over this period with monthly mean observed surface air 
temperatures (SATs) to produce La Nina-like ("ENSO C ool") and El Nino-like ("ENSO WA rm") conditions that are 
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A. ANN Ok ENSO warm SATfC) 0.15 

i i i : 



-3.6 -2.9 -2.3 -1.6 -.9 .5 .1 .3 .71.1 1.6 2 2.4 2.9 3.5 

B. ANN 0k AENSO cool SAT (°C) -0.10 



-3.6 -2.9 -2.3 -1.6 -.9-.5-.1 .3 .71.1 1.6 2 2.4 2.9 3.5 


Figure 1. Annual average surface air temperature anomalies (°C) used as boundary conditions for the (a) ENSO W arm and 
(b) ENSOcool simulations, relative to the preindustrial control mean. Average observed surface temperature anomalies 
[Reynolds and Smith, 1 994] associated with strong (x3) El Nino- (ENSOwarm) and La Nina-like (ENSOcool) conditions are 
determined from the NIN03.4 index and applied to SAT fields for the 0 k simulations as surface boundary conditions. Also 
shown are the locations of five defined regions of interest: the western (1), eastern (2) and central Pacific (3), eastern China 
(4), and northwest (NW) Pacific (5). Figure adapted from Lewis et at. [2013]. 

applied as surface boundary conditions for the preindustrial simulations. This technique produces surface 
temperature conditions that are broadly representative of positive and negative ENSO phases, although the 
regression was multiplied by a factor of 3 to ensure that the fingerprint of this type of climatic change is 
discernible throughout the hydrological cycle [Hoerling et at., 2001; Chung et at., 2014]. To put the 
temperature anomalies used into context, the corresponding NIN03.4 anomaly applied in our study is 
approximately equal to that occurring during the strong El Nino event during 1997-1998. Overall, we 
produce a preindustrial experiment suite with three ensemble members (0 k ENSO WAR m, 0 k ENSOcool. and 
0 k ENSO neutral ) that represent average, perpetual El Nino-, La Nina-, and neutral-like conditions. The surface 
temperature anomalies used for ENSO W arm and ENSO CO ol are shown in Figure 1. 

In addition, ENSO-like surface anomalies were also applied to the mid-Holocene for the various phases 
("6 k ENSOwarm" and "6 k ENSOcool") to investigate whether the fingerprint of ENSO-like variability on the 
hydrological cycle is sensitive to differing boundary conditions. For this experiment set, ENSO-like anomalies 
were applied as surface boundary conditions to mid-Holocene surface temperature and sea ice anomalies, 
which were calculated as differences between mean mid-Holocene and preindustrial conditions so that forced 
ENSO and forced mid-Holocene response could be separated. This results in a second experiment suite with 
three ensemble members — 6 k ENSO WA rm/ 6 k ENSOcool. and 6 k ENSO NEUTRA l. Although the amplitude and 
precise spatial pattern of ENSO activity in the mid-Holocene remains uncertain, we apply the same strong 
surface temperature forcing to both time slices to allow comparison of ENSO-like impacts on the tropical 
hydrological cycle under these different boundary conditions. These results are not presented for direct 
comparison with proxy records but rather as idealized experiments that help elucidate changing ENSO 
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A. Composite GPCP observed El Nino phase (DJF) 


C. Composite GPCP observed La Nina phase (DJF) 
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-20 -12 -4 


4 


12 


20 


Figure 2. Comparison of average December-February composite GPCP observational (a) El Nino precipitation anomalies 
with (b) simulated precipitation AENSOwarm anomalies (mm/d, relative to preindustrial simulation) and comparison of 
average December-February composite GPCP observational (c) La Nina precipitation anomalies with (d) simulated preci- 
pitation ENSOcool anomalies (mm/d, relative to preindustrial simulation). Observational El Nino (La Nina) precipitation 
anomalies (relative to a 1981-2005 climatology) are composites of 1982 and 1997 (1988 and 1999) (wherein sea surface 
temperature anomalies in the NIN03.4 region exceeded 2 standard deviations above (below) average for six consecutive 
months and represent strong El Nino (La Nina) events). 


teleconnection patterns and moisture source changes under varying boundary conditions. This approach also 
avoids relying on the model's ability to realistically capture intrinsic, freely evolving ENSO activity, which is 
inherently limited by model resolution and parameterizations, as well as our theoretical understanding of ENSO 
dynamics [Meehl et al., 2001]. The ENSO-like results should be interpreted in general terms, given the forcing 
applied corresponds to strong ENSO activity in a contemporary context. 

3. Evaluation of Model Performance 

In previous studies, isotope and precipitation distributions simulated with GISS ModelE-R have been 
systematically compared to observational records for the preindustrial [Schmidt et al., 2006, 2007] and 
evaluated against proxy records for the mid-Holocene [LeGrande and Schmidt, 2009]. Patterns of ModelE-R 
simulated 5 18 O p are highly correlated with observed fields, although rainfall in the western Warm Pool region 
is in excess of observed GPCP (Global Precipitation Climatology Project) patterns [ Schmidt et al., 2006]. 
Similarly, simulated mean Holocene East Asian 8 18 O p and surface ocean § 18 0 values correspond closely with 
those reported in regional proxy records, and the large-scale features of climatic change in the mid-Holocene 
are captured by the model [LeGrande and Schmidt, 2009]. 

In addition to the previous evaluation of the hydrological responses in these preindustrial and mid-Holocene 
model experiments, we evaluate the utility of the ENSO-like experiments for assessing changes in 
hydrological responses by comparing simulated spatial patterns with observed precipitation distributions 
(Figure 2). Simulated ENSO WARM and ENSO CO ol precipitation fields are compared to the GPCP gridded data 
sets [Adler et al., 2003], where the observed anomalies associated with strong El Nino conditions (i.e., in 
1982 and 1997) and strong La Nina conditions (i.e., in 1988 and 1999) were composited. The similarity of 
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observed (strong El Nino) and 
simulated (ENS0 WARM ) composite 
responses is seasonally variable; 
spatial anomalies are most similar in 
December-February (Figure 2) 
compared to during June-August 
(not shown). Similarly, observed and 
simulated ENS0 C ool phase anomalies 
are most alike in the boreal summer, 
although the simulated dry conditions 
over the central Pacific are notably zonal 
compared with observed. While the 
simulated precipitation responses in the 
ENSO-like experiment is generally 
stronger than those observed during 
these noteworthy ENS0 events, we 
consider these experiments to be 
insightful for examining potential 
changes in responses under varying 
boundary conditions. Flowever, as ENS0 
events typically reach their peak 
amplitude in the boreal winter 
(December-February (DJF)) [McGregor 
et al., 2013b], we consider modeled 
changes during this season only. 

In this study, model data are 
investigated in five regions (Figure 1, 
western, eastern and central Pacific, 
eastern China, and northwest (NW) 
Pacific) that are identified as sensitive 
to contemporary ENS0 variability, 
providing coverage of simulated 
hydrological change across the Pacific 
basin and encompassing the location 
of numerous proxy archives utilized in 
paleoclimate reconstructions of the 
Flolocene. These regions are not 
intended to provide exhaustive 
coverage of ENSO impacts, and hence, 
regions exhibiting significant ENSO- 
related climatological anomalies that 
do not encompass proxy archives, 
such as in the South Pacific and Indian 
Oceans, are not discussed explicitly. 

Simulated precipitation anomalies in 
these regions for ENS0 WARM and 
ENSOcool conditions are compared 
with observations in Table 2, showing 
that observed precipitation amount 
anomalies associated with strong 
observed El Nino and La Nina events 
are generally well captured in the 
model simulations in these regions. 
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Figure 3. Large-scale simulated (a, b) DJF precipitation (mm/d) and (c, d) 8 1s O p (%o) anomalies for the preindustrial 
simulation, showing ENSOwarm and ENSOcool changes relative to neutral conditions. All values reported are greater 
than 95% significant (Student's f test) given the year-to-year variability about the 5 year mean. 


Most notably, the simulated precipitation amount (3.6 mm/d) in the central Pacific is subdued relative to 
observed (9 mm/d), although only anomalies are considered in this study to account for consistent biases. 
Observed surface air temperatures and anomalies for ENSO WARM and ENSOcool conditions are also generally 
well captured in the model simulations. Simulated temperature anomalies are most disparate for the positive 
El Nino-like phase for the NW Pacific and eastern China regions, although it should be noted that these are 
compared to a limited number of observed years when other climate drivers may have also been regionally 
important. A detailed evaluation of simulated isotopic responses in these experiments is provided in the 
supporting information. 


4. Results 

The ENSO anomalies presented here represent the difference between the ENSO WARM or ENSOcool 
experiment and the ENSO neutral (ENSO neutral ) experiments. For example, using this notation, differences 
between simulated El Nino-like (La Nina-like) and neutral conditions in the preindustrial (Ok) experiment 
would be indicated as 0 k AENSO warm (0 k AENSOcool)- Similarly, we denote mid-Flolocene anomalies 
relative to preindustrial (0 k) values by A6k. Unless otherwise indicated, mid-Flolocene changes refer to mean 
ENSO neutral anomalies. Values reported here are greater than 95% significant (Student's t test), given 
year-to-year variability. 

4.1. ENSO-Like Hydrological Changes 

Simulated 8 18 O p and precipitation ENSO anomalies (AENSO warm and AENSOcool) for the preindustrial 
experiment set are shown in Figure 3 and summarized in Table 2. Under the warm sea surface temperature 
conditions in the tropical Pacific (characterizing the ENSO warm experiment), there is an increase in 
precipitation in the equatorial Pacific and decreases in the western Pacific around the Warm Pool. Simulated 
isotopic differences are spatially more variable than differences in precipitation amount between the 
ENSO warm and neutral experiments: ENSO-associated 5 18 O p anomalies extend to the high latitudes and 
display less spatial coherence than rainfall changes. Spatially variable 5 18 O p model responses, compared 
with precipitation changes, have been previously noted for other model experiments [ Schmidt et a!., 2005; 
Lewis et at., 201 0]. Generally, ENSO WARM § 18 O p is enriched (compared with neutral conditions) over the Warm 
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E. NW Pacific Ok ENSO neutral VSD (DJF) 



Figure 4. Simulated vapor source distributions (VSDs) for each analyzed region ((a) western Pacific, (b) central Pacific, 

(c) eastern Pacific, (d) eastern China, and (e) NW Pacific), showing the source for each region in the preindustrial simulation 
for the ENSO neutral phase. The solid rectangular boxes indicate the location of end-member precipitation. VSDs are 
unitless probability density function, normalized by the maximum probability density. 


LEWIS ETAL. 


©2014. American Geophysical Union. All Rights Reserved. 


7071 


^AGU Journal of Geophysical Research: Atmospheres 


10.1 002/201 3JD020961 


Pool region, extending through to the northern Indian Ocean. Conversely, over the eastern and central 
equatorial Pacific, there is a S 18 O p depletion during the ENSO WARM phase. 

There are areas over Australia, Africa, tropical South America, and the western USA in which distinct 5 18 O p 
responses to the ENSO WARM and ENS0 CO ol forcings are simulated, although no comparable differences in 
rainfall amount occur. Certain areas may be recording complex changes in hydrology, such as moisture 
source region, beyond simple deviations in rainfall amount alone [Vuilleand Werner, 2005; Schmidt et al., 2007; 
LeCrande and Schmidt, 2009]. Changes in moisture source regions to particular sites provide both an 
important control on 6 18 O p composition [Rozanski et al., 1 993] and a useful circulation diagnostic for 
considering regional-scale hydrological change [Lewis et al., 2010]. 

Regionally specific moisture source tracers demonstrate source changes during ENSO-like phases, relative to 
neutral conditions (Figures 4 and 5), and differences in regional hydrology occur across the Pacific basin 
between El Nino-like and La Nina-like events. Source region shifts associated with ENSO WARM conditions are 
generally stronger than the ENSO C ool, while the spatial patterns of change are distinct (Figure 5). That is, the 
ENSO warm anomalies, relative to neutral, do not necessarily mirror those associated with the ENSO C ool 
phase in all regions. An asymmetry in the ENSO-tropical rainfall teleconnection has been noted previously in 
response to a zonally symmetric ENSO forcing [Power et al., 2006; Chou and Lo, 2007]. The source region shifts 
associated with the ENSO CO ol experiment are largest in the central Pacific. The impact on the hydrological 
cycle of these ENSO-like surface temperature changes, even beyond the tropics, is also notable in the eastern 
China and northwest Pacific regions, although it should be reiterated that strong surface temperature 
anomalies were applied to these experiments as a forcing. In addition, the ENSO-phase source region 
changes are largely distinct from other climate modes, such as the seasonal cycle (not shown), although the 
amplitude of the moisture source response is comparable. 

4.2. Mid-Holocene Hydrological Changes 

In the mid-Holocene, simulated global DJF surface temperatures are 0.5°C cooler than the preindustrial 
(Annual (ANN), 0.35°C and June-August (JJA), 0.2 C). Climatological changes in the mid-Holocene, 
compared with the preindustrial, largely relate to insolation anomalies. In particular, due to changes in 
precession, there was reduced mid-Holocene insolation during DJF, resulting in cooler boreal winters and 
austral summers. In DJF, the spatial pattern of mean simulated temperature changes is characterized by 
cooling in the tropics, midlatitudes of the Northern Hemisphere, and southern high latitudes (Figure 6). 
Conversely, there is warming in the mid-Holocene in the high northern latitudes, compared with the 
preindustrial, due to increased wintertime insolation and the resulting, persistent reduction in sea 
ice [LeGrande and Schmidt, 2009]. 

Simulated mid-Holocene hydrological anomalies are less spatially coherent than surface air temperature 
changes. Precipitation amount anomalies are centered on the tropics, with mid-Holocene rainfall decreases 
simulated over the equatorial Pacific, southern Africa, South America, and Australia. Increases in precipitation 
are simulated across the Indian Ocean and Warm Pool region. The globally averaged mid-Holocene 
precipitation amount anomaly is 0.01 mm/d for DJF (ANN, —0.2 mm/d and JJA, —0.05 mm/d). The rainfall 
anomalies in the mid-Holocene, associated with orbital changes, are substantially smaller than those 
associated with the imposed ENSO-like forcings for the preindustrial (Figure 2), which are up to an order of 
magnitude higher. The simulated A6k 8 ls O p anomalies during DJF are more spatially widespread than the 
rainfall and surface air temperature changes, relative to the preindustrial. There are simulated isotopic 
enrichments throughout the high northern latitudes during the boreal winter and large-scale isotopic 
depletions across Eurasia and the Indian Ocean. The isotopic responses to insolation forcing in the 
mid-Holocene are similarly muted in comparison to the magnitude of the 5 18 O p anomalies during DJF 
associated with ENSO-like conditions (Figure 2). 

4.3. Mid-Holocene ENSO-Like Teleconnection Changes 

Comparison of S 18 O p and precipitation anomalies between the preindustrial and mid-Holocene indicates 
that the spatial patterns of change associated with ENSO-like phases are similar in each experiment (Figures 3 
and 7), with spatial fields highly correlated between the mid-Holocene and preindustrial (/?> 0.95). The 
ENSO-like precipitation amount and 8 ls O p changes are somewhat subdued in the mid-Holocene, in terms 
of both the magnitude of hydrological changes and their spatial extent. However, the precipitation source 
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Figure 5. Simulated vapor source distributions (VSDs) for each analyzed region ((a, f) western Pacific, (b, g) central Pacific, (c, h) 
eastern Pacific, (d, i) eastern China, and (e, j) NW Pacific), showing the difference (AENSO) between the ENSO W arm and neutral 
phases in the preindustrial control experiment (Figures 5a-5e). The AENSOcool differences in the preindustrial simulations, 
relative to neutral, are also shown (Figures 5f-5j). The solid rectangular boxes indicate the location of end-member precipitation. 
VSDs are unitless probability density function, normalized by the maximum probability density. 
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Figure 6. DJF seasonal average SAT (°Q, precipitation (mm/d), and 5 18 O p (%o) anomalies (A6k) for the mid-Flolocene, 
relative to the preindustrial simulation. All values reported are greater than 95% significant (Student's f test) given the 
year-to-year variability about the 5 year mean. Figure adapted from LeGmnde and Schmidt [2009], 


distributions for each of the analyzed regions (Figure 8) reveal a more complex relationship between 
ENSO-like expressions and the boundary conditions associated with the preindustrial and mid-Flolocene 
experiments. That is, source region difference between the mid-Flolocene and the preindustrial are less 
notable for the neutral ENSO conditions than when ENSO-like forcings are applied. The comparatively more 
pronounced impacts of ENSO CO ol and ENSO WARM , than neutral, conditions are a basin-wide phenomenon, 
impacting all regions analyzed here. 

For the western Pacific region, for example, there is a more westward (and localized) moisture source 
to the region of interest during the mid-Flolocene ENSO WA rm experiment than with equivalent surface 
temperature forcings applied during the preindustrial. Similarly, more localized moisture sources to 
the western Pacific region are simulated mid-Flolocene ENSO WARM experiment, with a decrease in 
rainfall sources from the equatorial Pacific, compared with the Warm Pool region. In addition, the spatial 
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Figure 7. As for Figure 3 but showing mid-Flolocene precipitation (mm/d) anomalies and 5 ls O p (%o) anomalies for (a, c) 
ENSO warm and (b, d) ENS0 CO ol changes, relative to neutral conditions. 


patterns of A6k anomalies for the ENSO C ool experiment do not necessarily mirror those for the 
ENSO warm experiment. For example, the differences between the mid-Holocene and preindustrial source 
distributions for the northwestern Pacific are weaker for the La Nina-like experiment suite than for the 
El Nino-like experiments. 

Comparing simulated surface wind changes between the mid-Holocene and preindustrial for the ENSO WARM 
and ENSO CO ol phases also reveals subtle but potentially important differences (Figure 9). Changes during 
the mid-Holocene include an outflow of surface air from eastern Asia, compared to the preindustrial, which 
shows landward airflow. The character of surface winds in the ENSO-like experiments is also notably different 
in the mid-Holocene compared to the preindustrial. For the ENSO WARM phases, particularly in the northern 
Pacific, flow is more zonally oriented during the preindustrial, and convergence along the equatorial Pacific is 
strong. Also, for the ENSO CO ol phase, the character of surface wind flow in the northwestern Pacific is distinct 
between the preindustrial and mid-Holocene, where convergence in the northern Pacific is weaker and 
shifted southwestward, which likely relates to modeled source region shifts (Figures 5 and 8). Hence, 
although the relationship between 5 1s O p and precipitation may be relatively stable for some areas under 
differing boundary conditions, the vapor source tracers and surface wind fields suggest that the interaction 
between ENSO-like variability and the large-scale ocean-atmosphere states (i.e., boundary conditions) may 
elicit nonlinear hydrological responses that warrant further investigation. 

5. Discussion 

5.1. Stability of Mid-Holocene ENSO-Like Expressions 

The assumption of stationarity in climatic relationships between regions underpins many paleoclimatic 
reconstructions. However, teleconnected climatic relationships are influenced by changes in the large-scale 
state of the coupled ocean-atmosphere systems that may have implications for paleoclimate research aiming 
to reconstruct changes in ENSO variability from information gleaned from one site alone [ Gallant et a/., 201 3]. As a 
simplistic assessment of potential differences in ENSO-like expressions between the preindustrial and 
mid-Holocene, we compared the relationship between 8 1s O p and precipitation amount across the different 
experiments (Figure 10). In certain regions, there is a clear distinction in the relationship between 5 18 O p and 
rainfall amount between the mid-Holocene and preindustrial for the various ENSO-like phases. For example, for 
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Figure 8. Simulated December-February precipitation source distributions for each region ((a) western Pacific, (b) central Pacific, (c) eastern Pacific, (d) eastern China, 
and (e) NW Pacific), showing the difference (A6 k) between the preindustrial control (0 kyr) and the mid-Flolocene (6 kyr) for each ENSO-like phase ((left) 
ENSO neutral , (middle) ENSO W arm- and (tight) ENSO CO ol)- The solid rectangular boxes indicate the location of end-member precipitation. VSDs are unitless 
probability density function, normalized by the maximum probability density. 


the northwestern Pacific region, the 8 1s O p -precipitation anomalies for ENSO WARM conditions (relative to neutral) 
are similar in both the preindustrial and the mid-Holocene (Figure lOe), although the ENSO C ool hydrological 
impact in this region is far more pronounced in the preindustrial and may relate to source region changes. 
Similarly, for the western and central Pacific region, the mid-Holocene, and preindustrial AENSO C ool anomalies are 
comparable, while the ENSO WARM phase anomalies are comparatively larger in the preindustrial than for 6 kyr. 
Next, for the eastern China region, the difference in the 8 18 O p -precipitaiton relationship between the mid- 
Holocene and preindustrial is comparable for the ENSO WARM and ENSO C ool phases (Figure lOd). 

Although this a simple metric of ENSO-like changes, it is clear that La Nina-like and El Nino-like hydrological 
expressions during the mid-Holocene will not necessarily follow those occurring during the preindustrial. In 
particular, the mid-Holocene hydrological (S 18 O p and precipitation amount) responses to ENSO WARM and 
ENSOcool phases are largely muted compared to the preindustrial, even though the same surface 
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Figure 9. Seasonal (DJF) surface wind speed anomalies (m/s) for the Pacific region for AENSOwarm and AENSOcool for 
the (a, c) preindustrial and (b, d) mid-Flolocene. Vectors indicate the direction of wind speed anomalies. 


temperature forcings are imposed for both time slices. While the simulated 8 1S 0 P and precipitation amount 
responses to ENSO-like forcings are subdued in the mid-Holocene, changes are simulated in other circulation 
diagnostics, such as the moisture source regions to the analyzed regions. When comparing source region changes 
during ENSO-like conditions during the preindustrial with those during the mid-Holocene (Figures 5 and 8), it is 
evident that ENSO-like hydrological expressions during the mid-Holocene are unique to that period. 

Furthermore, there are distinct patterns of hydrological changes associated with El Nino-like and La Nina-like 
conditions, and the application of ENSO C ool forcings does not produce climatological anomalies that 
represent impacts of equal magnitude but opposite sign of the ENSO WARM experiment. Rather, as previous 
research has also shown, the differing phases of ENSO are capable of producing nonlinear responses in 
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Figure 10. Comparison of 5 18 O p (%o) and precipitation amount (mm/d) ENSO-like phase anomalies for the preindustrial (0 k) 
and mid-Holocene (6 k) experiments for each defined region ((a) western Pacific, (b) central Pacific, (c) eastern Pacific, (d) 
eastern China, and (e) NW Pacific). The differing ENSO-like phases, relative to ENSO neutral , are designated by color (AENSO WA rm 
as red and AENSOcool as blue). Data points represent the areal average value in each of the regions, and arrows indicate the 
direction of statistically significant changes between the preindustrial and mid-Holocene for each phase. 
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precipitation, sea level pressure, and winds [Power et al., 2006; Conroy et al., 2013]. There are also distinct 
patterns of hydrological changes associated with the positive and negative ENSO phases between the 
preindustrial and mid-Holocene. Even in remotely teleconnected regions, there are changes associated with 
the ENSO-like phases that indicate nonlinear interactions between ENSO and the background state of the 
ocean-atmosphere system (i.e., the model boundary conditions). 

Why are expressions of ENSO-like variability distinct under different model boundary conditions, such as for 
the mid-Holocene? First, the complexity in responses to the imposed ENSO WARM and ENSO CO ol phases may 
relate partly to the differing mean climatic state that characterizes the mid-Holocene. The mid-Holocene 
tropical Pacific is generally more ENSO-like with larger temperature anomalies (relative to the preindustrial) 
occurring in the eastern Pacific rather than the western Pacific, despite a zonally symmetric insolation forcing 
(Figure 6). Hence, when surface temperature forcings are applied in distinct spatial regions (the central and 
eastern Pacific, compared to the western part of the basin), recognizably different responses may occur. 
Similarly, an increase in El Nino frequency has been simulated under projected future warming due to 
increased warming of the eastern equatorial Pacific, compared with the surrounding ocean [Power et al., 
2013; Caietal., 2014], Modern ENSO behavior is also sensitive to changes in wind stress that relate to surface 
temperature anomalies [Wittenberg, 2004], and changes in wind stress simulated here (Figure 9) may impact 
the hydrological cycle through convergence and ultimately convection. 

In addition, one of the predominant characteristics of ENSO is its connection to the seasonal cycle, with most 
recorded events reaching their maximum intensity in November-January. The phase-locking of ENSO 
variability to the annual cycle relates to the seasonal movement of the Intertropical Convergence Zone or 
zonal sea surface temperatures gradients, among other mechanisms [Ham et al., 2013]. Although the 
differences in boundary conditions between the mid-Holocene and the preindustrial are relatively small 
compared to other time slices, such as the Last Glacial Maximum, there are important perturbations in 
radiative forcings associated with this period [ Schmidt et al., 2004]. The largest changes are associated with 
varying insolation from the mid-Holocene to the present, primarily due to changes in precession, which 
results in a redistribution of seasonal insolation in the tropics. Specifically, tropical areas received more (less) 
irradiance during the second (first) half of the year and mid-Holocene Northern Hemisphere summers had 
substantially more incoming solar radiation than present (>25 W/m 2 over the season). The redistribution of 
seasonal insolation may impact tropical ocean-atmosphere dynamics in several ways, including through the 
perturbation of the zonal SST gradient and easterly trade winds [Clement et al., 1 999]. Also, changes in 
meridional SST gradients across the tropical Pacific, due to differential heating, may have impacted mid- 
Holocene upwelling and ENSO development [McGregor et al., 2013a]. 

5.2. Reconstructing Mid-Holocene ENSO Expressions 

The impact of the mid-Holocene orbital changes on the tropical Pacific and paleo-ENSO variability has not 
been established comprehensively from proxy records. Proxy ENSO-sensitive records indicate various 
possible changes in mid-Holocene ENSO characteristics in terms of the amplitude, variability, and spatial 
teleconnection patterns [ Donders et al., 2008]. Furthermore, some records suggest a relative cooling of 
eastern equatorial Pacific SST [ Koutavas et al., 2002], while others show relative warming [Lea et al., 2006]. In 
the western Pacific, S 18 0 records suggest similar or slightly cooler than modern SST, whereas concomitant Sr/Ca 
ratios indicate slightly warm SSTs [Tudhope et al., 2001]. Some fossil coral reconstructions indicate that 
there may have been reductions in ENSO variability in the mid-Holocene [Cobb etal., 2013], with integrated 
multiproxy reconstructions suggesting that recognizably "modern" El Nino activity beginning roughly 5 kyr 
[Donders et al., 2008]. 

Beyond the general constraints of reconstructing past climatic change from proxy archives that are often 
short or spatially limited, understanding paleo-ENSO variability is also complicated by the complex 
hydrological expressions of this large-scale phenomenon. Overall, these highly idealized experiments reveal 
inherent hydrological complexity associated with ENSO-like teleconnections that may impact interpretations 
of proxy ENSO-focused reconstructions. Our present model results, together with previous studies [Gallant et al., 
2013], highlight the importance of considering the potential nonstationarity of ENSO teleconnections when 
interpreting paleoclimate archives, especially those that rely exclusively on ENSO-related changes in precipitation. 
ENSO-like periodic patterns evident at proxy sites influenced by teleconnections may not necessarily reflect 
changes in ENSO-based frequencies but rather subtle changes in teleconnected climatological patterns 
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[Shulmeister et al., 2006]. In this case, interpretations of interannual-scale variability in climate proxies in terms of 
basin-wide ENSO expressions from one or two sites alone will be necessarily problematic. 

In addition, ENSO-like conditions impact many aspects of the climate system, including surface temperatures, 
wind speeds, and convective activity, as shown by the ENSO-forced circulation changes diagnosed by the 
modeled VSD tracers (Figure 8). In many regions, 5 18 O p variability integrates multiple climatic influences, such as 
changes in rainfall amount, temperature, and source regions, and hence, convoluted changes in ENSO-like 
expressions in the hydrological cycle are likely to make 5 18 0-climate interpretations inherently complex. For the 
long, continuous archives of climate change that are derived from § 18 O p variations, such as speleothems, those 
most useful for paleo-ENSO reconstructions may be those with the most direct climatic interpretations. That is, in 
locations like those situated in the equatorial Pacific, where simple, highly local S 18 0 climate interpretations of 
archives are likely to be most valid, the convoluting effects of ENSO-like changes, including moisture source 
regions, may be easier to constrain. Flowever, single site-based reconstructions should not be interpolated to 
regional or basin-wide scales for reconstructing past ENSO changes. Rather, to determine the spatial fingerprint of 
ENSO in the past, climatic responses must be synthesized across networks of reconstructions encompassing many 
sites. In addition, given the complexities of ENSO-phase hydrological expressions under differing boundary 
conditions, the interpretations of proxy records may benefit from the use of site-specific forward models that 
explicitly resolve the relevant, local hydrological response to particular ENSO conditions. 

5.3. Limitations of This Study 

The experiments utilized in this study are useful for revealing the complexity of hydrological responses to 
ENSO-associated surface temperature forcings. They do, however, include strong, idealized forcings, and we 
are limited in the conclusions that we can draw for locating proxy locations that are likely to record 
temporally stable ENSO impacts. Due to the difficulties in evaluating the utility of these experiments for real- 
world comparison, we considered only changes in the December-January season. Nonetheless, proxies can 
record signals in other seasons or can integrate annual averages. Ideally, a gridded proxy-based 
reconstruction of the tropical Pacific during the mid-Flolocene would benefit future model-based analyses of 
changing ENSO-like expressions. 

The comparison of the preindustrial and mid-Flolocene experiments was conducted to investigate potential 
changes under differing boundary conditions. We selected the mid-Flolocene as this is a well-constrained 
target for model-based studies [Schmidt et al., 2004], and we applied the same ENSO WARM and ENSO CO ol 
forcings as the preindustrial experiment suite. By applying the same forcings to both the preindustrial and 
mid-Flolocene experiments, we do not suggest that ENSO conditions during these two periods were 
indistinguishable. We also utilized a model with a relatively coarse horizontal resolution for the model experiments 
used in this study. Some model studies have demonstrated that precipitation regimes and other aspects of the 
hydrological cycle can be more realistically represented by models of high spatial resolution [Hirota et al., 2011]. 
Similarly, explicit model-based assistance in interpreting paleoclimate proxy records of ENSO-like variability would 
likely benefit for combining multiple isotope-enabled climate models [Conroy et al., 2013]. 

It is possible that we will only be able to understand long-term changes in ENSO dynamics with models with 
improved skill in producing physically realistic intrinsic, freely evolving ENSO variability. Flowever, problems 
have been identified in the capacity of GCMs to capture realistic cloud feedbacks, large-scale circulation 
patterns, and atmospheric waves that are necessary for elucidating ENSO expressions [Zhu et al., 2007; 
Dufresne and Bony, 2008]. Teleconnections in recent Coupled Model Intercomparison Project 5 (CMIP5) 
models show variable skill in different regions [Weare, 2013]. While models have generally improved 
compared to earlier versions, models remain most skillful in capturing observed precipitation and 
temperature anomalies in the tropics than in teleconnected regions. Furthermore, numerous GCMs do not 
simulate the phase-locking of ENSO to the annual cycle, due to erroneous zonal temperature gradients, 
coupling strengths, and thermocline depths [Ham et al., 2013]. 

6. Conclusions 

Using the water isotope- and vapor source distribution tracer-enabled GISS ModelE-R, we examined ENSO- 
like expressions in the hydrological cycle in various regions (eastern, western, central and northwestern 
Pacific, and eastern China) for multiple model time slices. Overall, these idealized model experiments reveal 
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hydrological complexity associated with ENSO-like conditions that may impact interpretations of proxy 
climate records. First, the regionally specific moisture source tracers incorporated into the model show that 
substantial changes in hydrology occur across the Pacific basin between the forced El Nino-like and La Nina- 
like events. The application of ENSO C ool forcings in the preindustrial simulation does not produce 
climatological anomalies that represent the equal but opposite impacts of the ENSO WARM experiment. 

Model results also demonstrate potential sensitivity in the interactions between ENSO-like phases and 
surface boundary conditions. Although the pattern of mid-Holocene changes for the differing ENSO phases 
are somewhat muted for rainfall amounts and 5 18 O p , changes are simulated in other circulation diagnostics in 
the mid-Holocene experiments, relative to the preindustrial. There are substantial changes in mid-Holocene 
moisture sources to the analyzed regions and surface wind fields, even in remote areas, which suggest 
complex interactions of the forced, idealized ENSO-like conditions with model boundary conditions. Overall, 
the mid-Holocene is characterized by a differing mean state compared to the preindustrial, with the largest 
difference in boundary conditions relating to orbital forcing. Differing mid-Holocene orbital configurations 
result in a redistribution of seasonal insolation in the tropics, which impacts surface temperatures, land-sea 
temperature contrasts, zonal circulation patterns, and wind stress. Although ENSO dynamics and 
teleconnections are known to be influenced by changes in the large-scale state of the coupled ocean- 
atmosphere systems and particularly the annual cycle, the assumption of stationarity in climatic relationships 
between regions and between ENSO phases underpins many paleoclimatic reconstructions. 

We suggest that the uncertainty surrounding the nonstationarity of ENSO teleconnections must be explicitly 
addressed in proxy-based paleo-ENSO reconstructions, particularly as model results demonstrate the potential 
sensitivity of ENSO-like phase impacts to surface boundary conditions. In particular, the simulations of the 
preindustrial and mid-Holocene demonstrate complex impacts on various aspects of the hydrological cycle, 
including moisture source regions, which may be integrated into proxy records (such as § 18 0 archives) and hence 
convolute reconstructions of paleoclimates. Ideally, our further understanding of paleo-ENSO changes will be 
guided by climatic responses synthesized across networks of reconstructions encompassing many sites, rather 
than from single sites alone. The synthesis of many sites will allow the spatial fingerprints of ENSO under varying 
boundary conditions (such as throughout the Holocene) to be determined more robustly and precisely. 

The scope of this study has been limited to investigating potential ENSO-like changes in the preindustrial and 
mid-Holocene time slices using idealized experiments and also been focused on isotopic variability. The 
experiments utilized here are useful for understanding hydrological responses to ENSO-associated temperature 
forcings under differing boundary conditions, but we reiterate that they do include strong, idealized forcings. As 
such, these simulations are limited in providing local-scale information for use in interpreting individual proxy 
records. Further model-based investigation of long-term ENSO dynamics, ideally using physically realistic freely 
evolving simulated ENSO variability, would be useful for targeting the location of paleoclimatic reconstructions 
toward regions with comparatively stable ENSO influences through time. In addition, the future use of multiple 
isotope-enabled climate models would help in interpreting paleoclimate proxy records of ENSO-like variability. 
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